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Highly-complex optical signal generation using electro-optical systems with 
non-linear, non-invertible transmission functions 
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We present a scheme whereby a static non-linear, non-invertible transmission function performed by the 
electro-optic Mach-Zehnder modulator produces highly complex optical chaos. The scheme allows the deter- 
ministic transformation of low-dimensional band-limited chaotic signals into much higher-dimensional struc- 
tures with broadband spectra and without using any delay elements or feedback. Standard benchmark tests 
show that all the considered complexity indices are highly increased due to this transformation in a con- 
trolled fashion. This mechanism allows the design of simple optoelectronic delayed oscillators with extremely 
complex chaotic output. 
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The generation of complex signals with broadband 
spectra is extremely important in diverse scientific and 
engineering applications. In particular chaotic systems 
have been implemented for this purpose using electri- 
cal and optical circuits. The construction of chaos-based 
encryption methods and the versatility of synchronizing 
chaotic systems have sparked a multitude of applications 
for electro-optical communication systems in the last two 
decades 1 . 

The well-known low-dimensional chaotic oscillators 
(Lorenz, Rossler, Chua, etc.) have been used for studying 
important dynamic properties like bifurcation, synchro- 
nization and spectral structures, but they lack perfor- 
mance when used for designing highly secure communica- 
tion systems. Therefore it has been necessary to use other 
schemes for producing more complex chaotic dynam- 
ics. Ikeda 2 introduced a system able to generate high- 
dimensional hyperchaos (multiple positive Lyapunov ex- 
ponents) that has been used as a model of electro-optical 
systems with delayed feedback 3 . These systems have 
shown a wide range of dynamics from bifurcations, multi- 
stability and chaotic breathers 4 in delayed optoelectronic 
oscillators (DOO) using the Mach-Zehnder modulator 
(MZM) as the source of nonlinearity. With these setups 
some advanced schemes of hyperchaos synchronization 
and optical communication systems 5 were developed. 

In this Letter we show a method to generate highly 
complex chaotic signals based on a robust theory. We 
show that using the MZM as a nonlinear static transfor- 
mation, we can strongly increase in a controlled way the 
dimension and other complexity measures of the output 
signals, including an enhancement of the bandwidth. We 
differentiate this scheme from others producing complex 
signals in that the source of the complexity is not given 



by any time delays, but by the topological properties of 
the transforming function. As the MZM is already an 
integral component of advanced chaotic communication 
schemes, our mechanism allows to improve the level of 
security and information-carrying capacity provided by 
common hyper chaotic systems. 

The theoretical foundation is based on the observation 
that some explicit deterministic functions produce ran- 
dom sequences of values 6 . These functions result from 
generalizing the solutions of standard chaotic mappings. 
The randomness here is defined in the sense of the im- 
possibility of writing future values of the sequences as 
single-valued relationships with past values (and vicev- 
ersa). For instance the function 



X n 



sm 2 (67rz n ), 



(1) 



where z and are real parameters, is the solution of the 
logistic map, X n+ i = 4X n (l — X n ), when z = 2. How- 
ever, it was demonstrated 7 that the sequences generated 
by this function cannot be written as a single- valued map 
of the form X n+i = f(X n , X n -i, X n - r +i), when z is 
not integer. Using this important finding we propose as 
model of our physical implementation the system 8 : 



X(t) = h(f(t)), 



(2) 



in analogy with the structure of Eq. (1). Deeper analyses 
have shown 7 ' 8 that the necessary conditions to produce 
complex dynamics in real physical systems using Eq. (2) 



Condition (i) : f(t) should be a non-periodic oscillat- 
ing function with repeating intervals of finite grow- 
ing behavior. 
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Condition (ii) : h(v) should be a non-invertible func- 
tion with many extrema such that h(v) = a pos- 
sesses several solutions, vi,V2, for v. 
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FIG. 1. Schematic diagram of the static non- linear, non- 
invertible system based on the MZM transmission function. 



For the physical implementation of this scheme we chose 
for f(t) a chaotic signal generated by a low-dimensional 
autonomous system. On the other hand, function h(y) 
will be implemented by the MZM, which intrinsically pro- 
duces the most ideal non-linear non-invertible function, 
a sinusoidal transformation. In Fig. 1 the experimen- 
tal setup shows a laser diode (1550 nm) coupled through 
mono-mode optical fiber into the MZM. The well-known 
structure of the MZM transmission is a non-linear func- 
tion of the applied voltage to the radio-frequency input 
(Vrf), (half- wave voltage V n = 5.2 V), and a DC voltage 
Vb to set the operating point on the transmission curve. 
Then the transmission can be written as: 



Po{Vrf) = 7^ cos 2 



ttVrf 



2V* 



(3) 



where P Q is the modulated output power, Pi is the power 
of the light entering the modulator, r] is a parameter ac- 
counting for the losses inside the device due to couplings 
and dissipations and <p is the bias term depending on 
Vb- The modulated light is coupled through another 
small section of mono-mode fiber into the photodetector 
(PIN diode) whose output voltage is sent to an oscillo- 
scope for recording. If in Eq. (3), we let Pi constant, 
set (j) = without losing generality, and let Vrf take 
a chaotic signal we obtain a topologically equivalent ex- 
pression to Eq. (2). Now, in order to ensure that the out- 
put will be of higher complexity than the input we have 
to fulfill Condition (ii). This means that the amplitude 
of the input has to span several lobes of the non-invertible 
transmission of Eq. (3). In the setup of Fig. 1 we use 
an amplifier of variable gain (Gi) that controls the mag- 
nitude of the input before entering into the MZM, and 
therefore we can control the number of spanned lobes of 
the nonlinear transmission. In the insert of Fig. (1), we 
amplified the input such that four lobes of the transmis- 
sion were spanned, or equivalent ly, a single output value, 
a, has a possibility of being produced by eight different 
input voltages i>i, v<i, ^8- This is the source of the in- 
determination and statistical independence of the output 
dynamics. This parameter G\ will give us another degree 
of freedom to study the relationship between the number 
of lobes of the transformation and the complexity of the 
output signal. This observation is key in the production 
of highly complex signals using the scheme of this Letter. 
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FIG. 2. Spectral characteristics of the chaotic signals, before 
(a) and (c), and after the transformation (b) and (d). For 
Duffing (upper) and Chua's (lower) oscillators. 
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FIG. 3. First return maps of the chaotic signals, before (a) 
and (c), and after the transformation (b) and (d). For 
Lorenz (upper) and Rossler (lower) oscillators. 



As a first application we use some well-known low- 
dimensional chaotic oscillators. The Chua, Rossler, 
Lorenz, and Duffing systems were implemented as elec- 
tronic circuits 9 . One component of each oscillator was 
amplified with the variable gain G\ (our complexity con- 
trol parameter) and inserted into the MZM. In Fig. (2), 
we show the spectral characteristic using the fast fourier 
transform (FFT) and first-return maps (FRM), compar- 
ing the input signals before and after the transformation. 
Interestingly the first observation is that the FFT of the 
transformed signals have a flattened noise-like spectra. 
The transformation increased their bandwidth, meaning 
that the original linear correlations in the chaotic signals 
are lost due to the transformation. This is in agreement 
with the theory from which this procedure is derived. 
If the signals lose correlation, the linear dependance be- 
tween future and past values is lost, and hence they ap- 
proach a random process in the sense defined above. This 
phenomenon is shown specifically for the chaotic dynam- 
ics of Duffing (Fig. (2-a,b)), and Chua's (Fig. (2-c,d)) 
circuits. We also show how the dimension, the ergodic- 
ity and uniform distribution of trajectories of the Rossler 
Fig. (3-c) and Lorenz Fig. (3-d) dynamics are increased 
in the two-dimensional projections given by the FRMs. 
The conclusions gathered from these figures can be but- 
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TABLE I. Complexity indices calculated from experimental 
time series before and after the transformation of the sys- 
tem of Fig. (1): False- neighbors dimension (FND), maximal 
Lyapunov exponent (A) and mutual information (MI). 
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tressed by a quantitative estimation of important param- 
eters that qualify the complexity of chaotic signals. There 
are several approaches to classify and study complexity 
of time series based on information theory and nonlinear 
analysis 10 . We use the very robust and professional free 
software package TISEAN 11 to obtain from the experi- 
mental time series the false nearest neighbor dimension, 
the maximal Lyapunov exponents and the mutual infor- 
mation. We chose these particular complexity indices 
because they are invariant measures and provide infor- 
mation about the underlaying dynamics rather than just 
simple statistics depending on the particular realizations 
used to calculate them 12 . The results of the calculations 
are compiled in Table (I). We can see that for every at- 
tractor transformed the dimension estimated by the false 
nearest neighbors (FND) method is strongly increased, 
obtaining the highest increment for the Rossler attrac- 
tor. However there was a weak increment in dimension 
for the Lorenz signal, which is due to the fact that the 
gain G\ used in this case allowed to span only two lobes 
of the transformation while for the Rossler signal there 
were six lobes. For Chua and Duffing the transformation 
used had four lobes. The maximal Lyapunov exponent, 
A, that measures the sensibility to initial conditions, or 
the chaoticity of the attractors are also significantly in- 
creased. Specially notable is the case where the Rossler 
attractor, the most quasi-harmonic oscillator, increases 
its maximal Lyapunov exponent 56 times when trans- 
formed. If the FFT figures above gave us an estimate 
of the elimination of linear correlations produced by our 
transformation, the mutual information (MI), accounts 
for non-linear correlations present in the time series 13 . 
When we calculate the MI on the transformed signals, 
the value r for the first zero in the MI is very close to 
zero, meaning this transformation eliminates almost ab- 
solutely any linear or non-linear correlations present in 
the original dynamics. 

We now turn to the question of the possible limitations 
of this method. That is, what happens when the origi- 
nal signal is already very complex and high-dimensional? 
Will the transformation produce the above results? To 
answer this we focus on delayed-feedback optoelectronic 
oscillators (DOO) which are known to generate some of 
the most complex chaotic behavior 14 . We implement as 
the source of the chaotic signal to be transformed a low- 
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FIG. 4. Schematic diagram of the static non-linear non- 
invert ible system based on the MZM-1 transmission function, 
transforming hyperchaotic signals from the DOO. 
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FIG. 5. First row: FFT of the hyperchaos generated by the 
DOO before (a) and after (b) the transformation. Second 
row: FRM of the hyperchaos before (c) and after (d) the 
transformation. 



frequency DOO as described by Murphy et al. 15 . Al- 
though this oscillator is very well described by many 
of the references of this Letter, we briefly mention our 
setup (see Fig. (4)): A 1550 nm laser diode (LD-2) 
with constant power (500 mW) lights through an op- 
tical fiber the MZM-2, the optical output of this device 
is converted to electrical voltage by a PIN photodetector 
(PD-2). Then this voltage is inserted into a digital signal 
processor board (DSP) TMDSDSK6713 of Texas Instru- 
ments. Inside the board we implement digital filtering to 
adjust the bandwidth of the chaotic signal in the range 
100 Hz-10 KHz as well as digital delay that we set in 2.6 
ms. This board is designed for audio frequencies with 
a maximum sampling rate of 96 KHz. After the signal 
is delayed and filtered by the DSP, it is converted back 
to analogue voltage and amplified by the constant gain 
G2 = 19.6. Finally this signal is inserted into Vrf of the 
MZM-2 (K = 4.5V) to close the loop. With the parame- 
ters set as above, the DOO performs hyperchaos with an 
FFT shown in Fig. (5-a). This signal is fed through the 



TABLE II. Complexity indices calculated from experimental 
time series before and after the transformation of the sys- 
tem of Fig. (4): False neighbors dimension (FND), maximal 
Lyapunov exponent (A) and mutual information (MI). 



FND 


A 


MI 


Attractor before after 


before after 


before after 


Hyperchaos 35 > 80 


1.5 3.0 


66.3 6.0 



TABLE III. KL dimension and entropy for different transfor- 
mations of the system of Fig. (4): KL dimension (KLD) and 
Shannon entropy (H). 
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amplifier with gain G\ = 2 into the transforming system. 
We would like to remark that the subsystems of Fig. (4) 
are two separate and independent entities with the upper 
system providing the hyerchaotic signal driving the lower 
system that performs the transformation. In the results 
shown in the first row of Fig. (5) we plot the FFT of the 
hyperchaos before (a) and after (b) the transformation. 
As in the low-dimensional case, we see that the correla- 
tions present in the original hyperchaos are lost in the 
transformation. In the second row we observe FRMs for 
the hyperchaos before (c) and after (d) the transforma- 
tion, the increase of dimension and ergodicity of orbits 
is again evident. Additionally and to compare with the 
low-dimensional study we obtain the same complexity 
indices shown in Table (II). We see that the dimen- 
sion, although originally high, is strongly increased up 
to the limits of calculations of the algorithms and soft- 
ware implementation given by the TISEAN. Also both 
the Lyapunov exponent and the mutual information are 
again significantly increased, meaning that transforming 
the hyperchaos produces a signal more hyperchaotic and 
totally uncorrelated. 

As a final test on the functionality of the method 
we added yet another rigorous analysis based on the 
Karhunen-Loeve (KL) decomposition 16 that has been 
implemented 17 on delayed electro-optical systems to as- 
sess how the time delays affect the complexity of the 
signals. This decomposition evaluates the dimension 
in terms of the number of principal orthogonal modes 
needed to reconstruct the signal and computes the Shan- 
non entropy from the eigenvalue spectra to quantify the 
dynamic complexity. In Table III we show the calcula- 
tions of the KL dimension and the Shannon entropy for 
different values of the gain (G\) of the dynamics gener- 
ated by system of Fig. (4). The first instance corresponds 
to the hyperchaos without transformation (G\ = 0). 
Then we increase by steps G\ to apply the transformation 
with increasing number of lobes from four to eight. The 



results show a monotone increase of the KL dimension 
and entropy, meaning that the complexity of the signal 
can indeed be controlled in this way. 

In summary we have shown a method to construct de- 
terministically very complex signals out of simpler in- 
puts, which can be readily implemented in the design of 
chaos-based optical communication systems. The signals 
produced have a wider bandwidth than those used to gen- 
erate them and have statistical properties of random pro- 
cesses with very high values of complexity quantifiers. In 
optoelectronic delayed oscillators, the chaotic complexity 
arises by virtue of delays in the system, in our scheme the 
complexity comes determined by the topological proper- 
ties of the transformation. Therefore an interesting ques- 
tion would be the relation arising between time delay and 
non-invertibility in the generation of very complex deter- 
ministic signals. We hope this work sheds some light and 
motivation in this direction. 
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